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Abstract 

We  report  experimental  observations  of  the  reaction  zone  structure  of  self-sustaining, 
cellular  detonations  propagating  near  the  Chapman-Jouguet  state  in  hydrogen- 
oxygen-argon/nitrogen  mixtures.  Two-dimensional  cross  sections  perpendicular  to 
the  propagation  direction  were  imaged  using  the  technique  of  planar  laser  induced 
fluorescence  (PLIF)  and  in  some  cases,  compared  to  simultaneously  acquired  schlie- 
ren  images.  Images  are  obtained  which  clearly  show  the  nature  of  the  disturbances 
in  an  intermediate  chemical  species  (OH)  created  by  the  variations  in  the  strength  of 
the  leading  shock  front  associated  with  the  transverse  wave  instability  of  a  propagat¬ 
ing  detonation.  The  images  are  compared  to  two-dimensional,  unsteady  simulations 
with  a  reduced  model  of  the  chemical  reaction  processes  in  the  hydrogen-oxygen- 
argon  system.  We  interpret  the  experimental  and  numerical  images  using  simple 
models  of  the  detonation  front  structure  based  on  the  “weak”  version  of  the  flow 
near  the  triple  point  or  intersection  of  three  shock  waves,  two  of  which  make  up 
the  shock  front  and  the  third  corresponding  to  the  wave  propagating  transversely 
to  the  front.  Both  the  unsteady  simulations  and  the  triple  point  calculations  are 
consistent  with  the  creation  of  keystone-shaped  regions  of  low  reactivity  behind  the 
incident  shock  near  the  end  of  the  oscillation  cycle  within  the  “cell” . 
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1  Introduction 


All  gaseous  detonation  waves  propagating  near  the  Chapman-Jouguet  velocity 
are  intrinsically  unstable.  The  instability  results  in  a  quasi-periodic  oscillation 
in  the  strength  of  the  leading  shock  wave  that  is  associated  with  a  set  of  weak 
shock  waves  propagating  transversely  to  the  main  front.  As  a  consequence, 
the  reaction  zone  region  is  spatially  nonuniform  and  unsteady.  The  three- 
dimensional  nature  of  the  detonation  front  was  discovered  by  White  [1]  with 
interferometry,  Denisov  and  Troshin  [2],  and  Shchelkin  and  Troshin  [3]  using 
soot  foils  and  has  since  been  studied  experimentally  by  many  researchers  using 
many  different  techniques. 

The  three-dimensional  complexity  of  the  front  creates  difficulties  for  flow  vi¬ 
sualization  techniques  that  integrate  through  the  flow  field.  To  overcome  this, 
many  detonation  structure  experiments  involving  flow  visualization  were  car¬ 
ried  out  in  narrow  channels:  facilities  in  which  at  least  one  dimension  is  on 
the  order  of  the  detonation  cell  width.  The  structure  of  the  triple  point  was 
studied  in  a  narrow  channel  using  schlieren  techniques  and  soot  foils  by  several 
researchers  [4-8] .  While  these  experiments  contribute  valuable  evidence  to  the 
study  of  detonation,  it  must  be  remembered  that  the  structure  of  a  detonation 
wave  is  altered  when  it  propagates  in  a  narrow  channel,  as  shown  by  Strehlow 
and  Crooker  [7]. 

Schlieren  images  of  detonations  in  non-marginal  cases  have  also  been  obtained 
by  Takai  et  al.  [9],  among  others.  In  the  case  of  mixtures  with  regular  structure, 
schlieren  images  show  similar  shock  configurations  to  those  observed  in  narrow 
channels,  when  the  transverse  waves  were  orthogonal  with  respect  to  the  side 
walls  and  the  optical  path.  In  nitrogen-diluted  mixtures,  the  schlieren  images 
also  show  a  more  irregular  cell  structure  as  observed  on  soot  foils. 

Direct  experimental  visualizations  of  chemical  species  in  a  detonation  wave 
are  more  limited.  Dormal  et  al.  [10]  combined  OH  emission  spectroscopy  with 
pressure  traces  and  soot  foils  to  measure  the  induction  time  through  a  cell. 
Andresen  et  al.  [11]  used  the  method  of  planar  laser  induced  predissociated 
fluorescence  (LIPF)  in  an  effort  to  visualize  a  cross  section  of  the  OH  con¬ 
centration  behind  a  detonation  front  in  a  hydrogen-oxygen-argon  mixture. 
Although  OH  images  were  obtained,  there  was  not  sufficient  spatial  resolution 
to  obtain  any  information  other  than  the  average  spatial  period  of  the  oscil¬ 
lations  in  the  OH  emission  at  the  detonation  front.  One  major  difficulty  that 
this  study  identified  was  the  substantial  amount  of  natural  fluorescence  that 
occurs  in  detonations  and  the  importance  of  reducing  this  in  the  final  image  in 
order  to  obtain  high  signal-to-noise  ratios.  Rayleigh  scattering  from  a  planar 
laser  sheet  was  used  by  Anderson  and  Dabora  [12]  to  image  a  cross  section 
of  the  density  field  behind  detonations  propagating  in  hydrogen-oxygen-argon 
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mixtures.  Shock  front  shapes  and  quasi-periodic  spatial  fluctuations  in  den¬ 
sity  were  observed  that  correlated  with  previous  numerical  and  experimental 
studies.  These  results  confirmed  the  density  variations  inferred  from  previous 
experimental  investigations  and  predicted  by  simulation.  However,  there  was 
no  insight  provided  into  the  chemical  processes  taking  place  in  the  reaction 
zone. 

More  recently,  LIPF  was  applied  with  greater  success  to  fast  deflagrations  and 
detonations  in  hydrogen-air  mixtures  by  Gerlach  [13]  and  Eder  [14].  Images  of 
the  chemical  reaction  front  with  sufficient  spatial  resolution  and  signal-to-noise 
ratio  were  achieved  by  using  a  tunable  KrF  laser  with  a  pulse  energy  of  450  mJ 
exciting  the  OH-radical  at  a  wavelength  of  248  nm.  A  set  of  filters  transmits 
the  fluorescence  signal  between  wavelengths  of  288  and  303  nm.  In  a  facility 
with  a  diameter,  D,  of  66  mm,  Eder  examined  a  variety  of  fast  combustion 
processes  in  hydrogen- air  mixtures  with  an  equivalence  ratio  between  0.4  and 
2.4  at  atmospheric  pressure.  For  marginal  detonations  with  X/D  ~  1.6  where 
A  is  the  detonation  cell  width,  an  irregular  reaction  front  containing  some 
unreacted  regions  was  observed  behind  the  leading  shock  wave.  In  detonations 
with  X/D  <  1,  the  reaction  front  looked  rather  smooth  and  it  was  difficult  to 
observe  evidence  of  the  cellular  structure. 

The  application  of  LIPF  and  Rayleigh  scattering  appeared  to  be  limited  for 
studies  of  fully  developed  detonations  given  the  relatively  poor  results  in  com¬ 
parison  to  the  success  of  LIF  methods  in  elucidating  the  reaction  zone  struc¬ 
ture  of  turbulent  flames  [15]  and  shock- induced  combustion  [16-18].  LIPF  is  a 
relatively  low  yield  process  and  it  was  apparent  that  the  signal-to-noise  ratio 
was  a  significant  issue  for  this  application.  Rayleigh  scattering  is  an  extremely 
low  yield  method  and  also  very  difficult  to  use  in  most  detonation  tubes  be¬ 
cause  of  the  large  numbers  of  particulates  and  soot  buildup  due  to  use  with 
hydrocarbons  or  soot  foil  diagnostics.  A  more  promising  technique  for  investi¬ 
gating  fully  developed  detonations  is  nonresonant  LIF,  pumping  and  detecting 
at  different  frequencies.  By  carefully  choosing  the  parameters  of  the  detona¬ 
tion  and  LIF  system,  we  felt  that  it  would  be  possible  to  overcome  previous 
difficulties  and  image  the  details  of  the  reaction  zone  structure. 

In  the  present  study,  we  build  on  the  extensive  prior  work  on  detonation 
structure  and  apply  both  classical  techniques  of  soot  foils  and  schlieren  imag¬ 
ing  alongside  the  newer  method  of  Planar  LIF  or  PLIF  to  visualize  the  OH 
concentration.  We  are  concerned  with  cases  of  multi-headed  waves  that  have 
an  instability  structure  that  is  fully  three-dimensional  and  perform  experi¬ 
ments  in  a  test  section  with  a  large  (150  mm  on  a  side)  square  cross  section. 
In  this  paper,  we  report  on  experiments  with  mixtures  of  stoichiometric  H2- 
02  diluted  with  either  Ar  or  N2.  After  describing  our  experimental  setup,  we 
give  the  experimental  results,  then  describe  numerical  simulations,  and  com¬ 
pare  our  experiments  to  the  simulations.  We  carry  out  a  detailed  analysis  of 
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the  triple  point  configuration  and  consider  the  implications  of  our  study  to 
the  broader  issues  connected  with  the  role  of  transverse  waves  in  detonation 
propagation. 


2  Experimental  Setup 


The  experiments  were  carried  out  using  a  7.3  m  long,  280  mm  diameter  deto¬ 
nation  tube  [19]  attached  with  a  “cookie-cutter”  to  the  150  mm  square  test  sec¬ 
tion  described  in  Kaneshige  [20].  The  facility  is  shown  in  Fig.  la.  Detonations 
are  initiated  with  an  exploding  wire  and  a  short  section  of  acetylene-oxygen 
mixture  injected  immediately  prior  to  the  wire  explosion.  The  acetylene-oxygen 
driver  is  injected  through  a  distributed  manifold  surrounding  the  wire.  The 
initiation  system  results  in  direct  initiation  with  a  highly  repeatable  opera¬ 
tion  and  detonation  velocities  within  0.5%  of  the  Chapman- Jouguet  value. 
Detonation  pressure  and  velocity  is  measured  by  a  series  of  pressure  trans¬ 
ducers  located  along  the  length  of  the  tube  and  in  the  test  section.  The  flow 
may  be  visualized  through  windows  located  approximately  7.5  m  from  the 
initiation  point.  Stoichiometric  mixtures  of  hydrogen  and  oxygen  diluted  with 
argon  or  nitrogen,  initial  pressure  of  20  kPa,  and  nominal  initial  temperature 
of  25° C  have  been  examined  in  the  present  study.  Schlicren  images,  pressure 
histories,  and  soot  foil  records  indicate  the  detonations  are  fully  developed, 
self-sustaining  waves  within  the  test  section. 

The  detonation  reaction  zone  has  been  visualized  by  exciting  OH  fluorescence 
at  about  284  nm  with  a  planar  light  sheet  and  imaging  the  emitted  light 
with  an  intensified  charged-coupled  device  camera.  The  development  of  the 
PLIF  visualization  method  and  the  initial  experimental  results  are  described 
in  detail  by  Pintgen  [21],  Simultaneous  observations  with  a  laser  schlieren 
system  have  also  been  used  in  some  tests. 

The  UV  light  for  excitation  of  the  OH  molecules  is  produced  by  frequency  dou¬ 
bling  the  output  of  an  excimer-pumped  dye  laser.  The  excimer  laser  (Lambda 
Physik,  Model  Cornpex  102,  XeCl  filling)  produces  a  17  ns  pulse  of  broad¬ 
band  light  centered  at  308  nm  with  a  total  energy  of  about  300  mJ.  The  dye 
laser  (Lambda  Physik,  Model  Scanmate  2E)  uses  Coumarin  153  dye  and  a  fre¬ 
quency  doubler  to  produce  narrowband  UV  light  with  a  total  energy  of  about 
6.5  mJ.  The  frequency  of  the  dye  laser  was  tuned  to  284.008  nm,  which  is 
located  half-way  between  two  OH  transitions:  a)  A2E+  <— X2n(l,0)  Q2(8)  at 
284.009  nm;  b)  A2X+  ^X2n(l,0)  Qi(9)  at  284.007  nm.  The  choice  of  the  two 
transitions  was  based  on  computations  with  the  LIFBASE  program  [22]  and 
repeated  experimental  trials.  The  offset  of  the  dye  laser  frequency  calibration 
was  determined  by  scanning  in  frequency  while  measuring  the  fluorescence 
yield  from  a  hydrogen  diffusion  flame  in  air,  comparing  the  results  to  sirnu- 
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lated  scans  with  LIFBASE.  The  output  of  the  frequency  doubler  was  formed 
into  a  sheet  by  using  a  combination  of  cylindrical  (focal  length  -25.4  mm)  and 
spherical  lenses  (focal  length  1000  mm).  The  final  usable  height  was  about 
80  mm,  and  the  thickness  of  about  0.3  mm  was  determined  by  burns  on  ther¬ 
mopaper.  The  light  sheet  entered  the  test  section  through  a  quartz  window  at 
the  end  wall  of  the  test  section. 

The  induced  fluorescence  emerged  through  a  quartz  window  at  90  degrees  to 
the  light  sheet  (Fig.  lb)  and  was  imaged  on  to  a  slow-scan  576x384  pixel 
ICCD  (Princeton  Instruments  ITE/ICCD-576)  camera  through  a  bandpass 
filter  with  a  centerline  of  313  nm  and  10  nm  FWHM.  The  detected  radiation 
is  expected  to  be  primarily  due  to  (1,1)  transitions  [23]  near  315  nm.  Detection 
in  this  spectral  region  minimizes  problems  such  as  fluorescence  trapping  and 
ensures  that  the  images  will  not  be  contaminated  by  Mie  scattering  or  stray 
reflections  from  the  test  section  of  the  incident  light  sheet.  An  image  was 
formed  by  a  105  mm  f/4.5  Nikkor  UV-transmitting  lens  on  to  the  microchannel 
plate,  which  was  gated  with  a  100  ns  pulse  of  800  V.  The  gate  signal  to 
the  camera  was  synchronized  with  the  laser  firing  pulse,  and  the  laser  firing 
was  coordinated  with  the  detonation  front  location  by  using  a  programmable 
delay  generator  triggered  from  a  pressure  transducer  located  upstream  of  the 
visualization  section.  The  signal  from  the  pressure  transducer  is  processed 
through  a  latching  edge-detection  trigger  circuit  to  eliminate  false  triggers 
from  the  capacitor  discharge  system  used  to  initiate  the  detonation. 

Careful  adjustment  of  the  timing  is  essential  to  obtaining  OH  images  since  the 
excited  state  population  quickly  decays  and  the  camera  gate  time  must  also  be 
limited  to  maintain  a  high  signal-to-noise  ratio  in  the  image.  There  are  very 
high  concentrations  of  OH  behind  a  detonation  front,  and  natural  fluorescence 
is  a  significant  source  of  noise.  Due  to  the  short  life  time  of  the  excited  state, 
the  laser  pulse  must  occur  within  the  camera  gate  time  in  order  to  obtain  an 
image  of  the  laser  induced  fluorescence.  In  the  present  case,  the  operation  of 
the  pulse-forming  circuit  and  thyratron  switch  in  the  excimer  laser  resulted 
in  a  delay  of  about  1200  ns  between  the  trigger  pulse  and  the  light  output.  A 
jitter  of  up  to  200  ns  can  occur  in  the  function  of  the  thyratron,  requiring  the 
timing  to  be  tested  and  adjusted  prior  to  each  detonation  experiment.  The 
camera  gate  pulse  was  set  to  occur  40  ns  prior  to  the  assumed  nominal  laser 
function  time. 

For  simultaneous  schlieren  and  PLIF  images,  the  Q-switched,  flash-lamp  pumped 
ruby  laser  described  by  Akbar  [19]  was  used  to  create  a  short  (<  50  ns)  light 
pulse  at  693  nm.  The  light  was  formed  into  a  collimated  beam  about  150  mm 
diameter,  passed  through  the  test  section,  and  then  imaged  by  a  combination 
of  mirrors  and  lenses  onto  high-speed  instant  (Polaroid  3000  ISO)  him.  The 
time  between  acquiring  the  schlieren  and  PLIF  images  is  measured  by  pho¬ 
todiodes  to  be  less  than  600  ns.  The  ICCD  camera  is  placed  just  outside  the 
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optical  path  of  the  schlieren  system  and,  therefore,  views  the  laser  sheet  not 
normally  but  at  an  angle  of  12  °(Fig.  1),  causing  the  PL1F  image  to  be  slightly 
distorted.  Images  of  a  test  pattern  and  postprocessing  are  used  to  correct  this 
distortion  and  also  to  align  the  PL1F  and  schlieren  image. 


3  Experimental  Results 


One  OH  image  is  obtained  for  each  experiment.  A  selection  of  these  images  is 
shown  in  Figs.  2-4;  the  image  is  identified  by  the  test  number.  The  height  of 
each  image  is  approximately  80  mm  and  the  detonation  wave  is  traveling  from 
left  to  right  in  all  images.  The  detonation  front  instability  exhibits  seemingly 
random  variations  in  the  phasing  of  the  transverse  waves.  For  this  reason,  it 
is  not  possible  to  use  repeated  experiments  to  construct  a  coherent  time  se¬ 
quence.  It  is  also  important  to  note  that  these  images  are  qualitative  rather 
than  quantitative  since  we  have  not  corrected  for  the  effect  of  the  laser  beam 
nonuniformity  or  molecular  effects  such  as  quenching  or  time-dependence  of 
the  excited  state  population.  The  main  value  of  the  present  results  is  in  re¬ 
vealing  the  geometry  of  the  OH  fronts  and  enabling  an  interpretation  in  terms 
of  gas  dynamics  and  chemical  reactions. 

A  distinct  OH  concentration  front  is  visible  in  all  cases.  Simulations  of  the 
reaction  zone  behind  the  shock  front,  discussed  subsequently,  indicate  that 
the  origin  of  the  front  is  clue  to  the  sharp  rise  in  the  OH  concentration  at 
the  end  of  the  induction  zone  or  region  of  radical  chain  reaction  behind  the 
main  shock  front.  The  fluorescence  intensity  decreases  rapidly  with  increasing 
distance  downstream  (to  the  left  in  the  images)  from  the  front  and  there  is 
significant  radiation  detected  from  only  the  first  20-30  mm  behind  the  front. 
This  can  probably  be  explained  by  the  dependence  of  the  PLIF  signal  inten¬ 
sity  on  the  OH  number  density  profile  through  the  reaction  zone.  The  OH 
number  density  rises  sharply,  reaches  a  maximum,  and  then  decreases  with 
increasing  distance  from  the  shock,  as  shown  by  a  ZND  calculated  profile  in 
Fig.  5.  A  temperature  profile  is  also  shown.  Other  parameters  influence  the 
spatial  variation  of  the  collected  signal  magnitude:  the  transition  probability 
of  absorbance  and  spontaneous  emission,  variations  in  the  X2n  initial  state 
population  with  the  local  thermodynamic  state,  and  quenching  of  the  excited 
A2£+  state  by  water  [24],  The  absorption  of  the  light  sheet  energy  by  OH  can 
also  gradually  decrease  the  fluorescence  intensity  behind  the  shock.  Due  to  the 
absence  of  large  shadows  and  caustics,  refraction  of  the  incident  light  associ¬ 
ated  with  the  large  refractive  index  variations  near  the  absorption  lines  seems 
to  have  a  minor  influence  on  the  PLIF  signal.  The  efficiency  of  the  collection 
optics  will,  of  course,  influence  the  overall  PLIF  signal  intensity  [25]. 

The  geometry  of  the  OH  front  is  striking  and  some  common  features  stand  out. 
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A  single  and  particularly  pronounced  feature  with  a  keystone  shape  is  seen  in 
shot  1419.  Similar  features  appear  in  other  images.  An  important  point  to  note 
while  examining  these  images  is  that  the  light  sheet  may  not  have  intersected 
the  cellular  structure  orthogonally,  as  can  be  seen  from  soot  foils  taken  at 
the  test  section  end  wall  (Fig.  6).  That  is,  although  the  light  sheet  is  aligned 
perpendicularly  to  the  test  section  walls,  the  detonation  cellular  structure 
may  not  be  orthogonal  to  these  walls,  potentially  resulting  in  some  three- 
dimensional  effects  in  the  images.  In  the  facility  used  in  these  experiments, 
the  detonation  transitions  from  a  tube  with  a  round  cross  section  to  a  square 
test  section  via  a  cookie-cutter.  This  change  in  geometry  may  influence  the 
orientation  of  the  cell  structure.  Fig.  7  shows  schematically  how  such  non- 
orthogonal  intersections  may  result  in  features  similar  to  those  observed  in 
some  PLIF  images. 


In  Figs.  2-4,  the  shock  location  is  not  known  but  it  can  be  inferred  from  the 
simulations  and  simultaneous  schlieren  and  PLIF  images  to  be  just  to  the  right 
of  the  OH  emission  front.  We  anticipate  that  the  distance  from  the  shock  to 
the  OH  front  will  be  a  strong  function  of  shock  strength  and  the  clearly  visible 
sudden  changes  or  jumps  in  the  location  of  the  front  must  be  associated  with 
the  instability  of  the  detonation  wave. 


This  is  confirmed  by  the  simultaneous  schlieren  and  OH  PLIF  images  shown 
in  Fig.  8.  Despite  the  masking  effect  of  the  three-dimensional  nature  of  the 
transverse  waves,  it  is  clear  that  the  changes  in  the  OH  concentration  front 
location  can  be  correlated  directly  with  the  density  held  changes  visualized  by 
the  schlieren  images.  The  variation  in  the  OH  concentration  front  location  with 
shock  velocity  is  confirmed  by  numerical  simulations  described  subsequently. 


The  mixtures  chosen  in  these  experiments  were  intended  to  investigate  both 
regular  (Ar-diluted  mixtures)  and  more  irregular  (N2-diluted  mixtures)  cell 
structure  to  determine  if  there  were  corresponding  differences  in  the  chemical 
reaction  zone  structure  and  particularly,  if  differences  could  be  seen  in  the 
mechanism  by  which  combustion  occurs.  Previous  work  [26]  has  suggested 
that  regular  mixtures  react  by  auto-ignition  following  shock  compression,  while 
the  reaction  zone  of  irregular  mixtures  is  affected  and  possibly  sustained  by 
turbulence-producing  shock  wave  interactions  with  other  shock  waves  or  shear 
layers.  We  hope  to  be  able  to  investigate  these  issues  by  directly  visualizing 
the  reaction  zone  structure  of  mixtures  of  varying  cellular  regularity.  A  list  of 
mixture  parameters  is  given  in  Table  1.  Sooted  foils  from  a  280  mm  diameter 
detonation  tube  were  used  to  measure  the  cell  width  in  these  mixtures  and 
illustrate  the  differences  in  cellular  regularity  of  the  mixtures  chosen  (Fig.  9). 
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4  Numerical  Simulation 


Numerical  simulations  of  two-dimensional  detonations  traveling  in  a  6.0  cm 
channel  at  the  Chapman-Jouguet  velocity  in  2H2-02-7Ar  mixtures  at  an  ini¬ 
tial  pressure  of  6.7  kPa  have  been  carried  out  by  Eckett  [30].  These  simula¬ 
tions  used  a  reaction  model  and  thermochemistry  based  on  a  detailed  reaction 
mechanism.  The  chemical  reactions  were  represented  by  an  approximate  four¬ 
dimensional  manifold  using  the  ILDM  technique  augmented  with  an  induction 
time  manifold  for  the  initiation  stage;  details  are  provided  in  Eckett  [30].  The 
numerical  simulations  used  adaptive  mesh  refinement  and  an  approximate 
Riemann  solver  (Roe-Glaister)  to  provide  high  resolution  and  accuracy.  The 
ILDM  reduced  mechanism  and  the  numerical  solution  technique  were  vali¬ 
dated  against  several  test  cases  and  previously  published  results  using  the  full 
detailed  chemistry. 

Simulated  schlieren  and  OH  mass  fraction  images  are  shown  in  Fig.  10.  A 
close-np  of  the  OH  mass  fraction  and  the  key  features  of  the  detonation  front 
are  shown  in  Fig.  11.  Several  features  can  be  observed  in  the  close-up  of  the  OH 
mass  fraction.  First,  there  are  the  obvious  OH  fronts  closely  following  curved 
portions  (Mach  stems)  of  the  leading  shock  wave.  Second,  on  the  centerline 
between  the  two  Mach  stems,  there  is  a  keystone-like  feature  which  is  bounded 
by  the  lead  shock,  the  shear  layer,  and  the  OH  front  behind  the  incident  wave 
(the  geometrical  configuration  is  discussed  in  more  detail  in  the  following 
section).  Third,  the  decrease  in  mass  fraction  along  the  downstream  edge  of 
slip  lines  can  also  be  seen.  There  is  a  striking  qualitative  similarity  between 
the  overall  appearance  of  the  numerical  simulation  of  the  OH  mass  fraction 
and  the  experimental  images  presented  in  Fig.  3.  Keystone  features  are  present 
in  Fig.  2  but  since  the  reaction  zone  is  much  shorter,  they  are  not  as  distinct 
as  in  the  85%  argon  case.  The  features  in  Fig.  4  are  much  more  irregular 
and  pronounced  than  in  the  argon-diluted  cases.  This  is  consistent  with  the 
well-known  observations  of  the  irregularity  of  the  cellular  structure  in  nitrogen- 
diluted  cases  as  compared  to  the  regular  structure  with  large  amounts  of  argon 
dilution. 

It  is  important  to  note  that  the  three-dimensional  nature  of  the  experimen¬ 
tal  studies,  the  nonuniformity  of  the  light  sheet,  nonideal  effects  in  the  flu¬ 
orescence  process,  and  response  of  the  camera  will  be  important  factors  in 
determining  the  appearance  of  the  OH  front  in  the  image.  Given  all  of  these 
issues,  it  is  not  surprising  that  the  simulations  do  not  exactly  reproduce  the 
observed  shapes.  Nevertheless,  the  results  are  similar  and  provide  direct  evi¬ 
dence  of  how  the  chemical  reaction  processes  are  affected  by  the  detonation 
front  instability. 

Understanding  the  keystone  regions  requires  investigating  in  detail  the  struc- 


ture  of  the  reaction  zones  behind  the  leading  and  transverse  shock  wave  sys¬ 
tems  that  make  np  the  detonation  front.  Detailed  analysis  of  the  simulations 
can  be  used  to  obtain  the  variation  of  the  leading  shock  front  pressure  and  ve¬ 
locity  through  the  cell;  the  results  are  shown  in  Fig.  12.  The  pressure  and  shock 
velocity  variation  indicate  that  the  leading  shock  front  velocity  is  a  maximum 
value  of  about  1.3  Ucj  at  the  beginning  of  the  cell  (x  =  xQ ),  dropping  sharply 
to  about  0.93  Ucj  by  the  middle  of  the  cell  ( x  —  xQ=  0.5 L)  and  remaining 
nearly  constant  at  this  value  until  the  collision  of  the  transverse  waves  causes 
a  rapid  increase  in  the  shock  front  velocity  at  the  end  of  the  cell  ( x-xQ  =  L). 
The  observed  range  of  leading  shock  velocities,  0.91  <  U /Ucj  <  1.3,  is  similar 
to  that  computed  by  Oran  et  al.  [31]  but  smaller  than  has  been  experimen¬ 
tally  observed  by  previous  researchers  [4,32,7]  who  studied  marginal  waves,  for 
which  the  excursions  in  lead  shock  velocities  are  apparently  much  stronger.  In 
addition,  Gamezo  et  al.  [33]  calculate  that  the  amplitude  of  oscillation  of  the 
leading  shock  velocity  increases  with  increasing  mixture  activation  energy. 

One  key  issue  is  the  variation  of  the  rate  of  reaction  behind  the  shock  front.  A 
common  approach  to  this  issue  is  to  use  a  quasi-steady  estimate  of  the  reaction 
zone  structure  based  on  the  instantaneous  value  of  the  leading  shock  velocity. 
Shown  in  Fig.  13  is  an  estimate  of  the  induction  time  (based  on  the  location 
of  the  maximum  energy  release  rate)  for  the  simulations  shown  in  Fig.  10. 
The  quasi-steady  approach  ignores  the  time-dependence  or  unsteadiness  of  the 
shock  front  and  the  possible  effect  on  the  reaction  zone  structure.  The  role  of 
unsteadiness  in  detonation  initiation  has  been  considered  by  Eckett  et  al.  [34], 
They  show  that  the  effect  of  a  decaying  wave  is  to  increase  the  reaction  zone 
thickness,  and  the  amount  of  increase  is  a  strong  function  of  the  rate  of  decay 
of  the  leading  shock  front.  If  the  shock  front  decay  time  is  less  than  a  critical 
value,  then  Eckett  et  al.  [34]  find  that  the  unsteadiness  can  actually  quench  the 
reaction  and  cause  the  energy  release  to  become  completely  decoupled  from  the 
shock  front.  Previous  experimental  and  numerical  studies  on  detonations  have 
suggested  that  this  occurs  in  some  mixtures,  resulting  in  unreacted  “pockets” 
behind  the  detonation  front. 


5  Unreacted  Gas  Pockets 


An  unburnt  gas  pocket  is  defined  here  as  a  region  of  unreacted  gas  that  be¬ 
comes  detached  from  the  main  detonation  front  after  the  collision  of  two  triple 
points.  It  should  not  be  confused  with  the  unreacted  gas  in  the  induction  zone 
behind  the  incident  wave  which  results  in  the  keystone  features  observed  in  the 
PLIF  images.  Unburnt  pockets  are  speculated  to  persist  to  a  distance  on  the 
order  of  a  cell  dimension  downstream  of  the  main  front  and  to  be  resolvable 
at  the  scale  of  the  cell  width. 
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“Islands”  of  unreacted  gas  that  become  isolated  downstream  of  the  main  det¬ 
onation  front  were  first  discussed  by  Subbotin  [8]  in  his  schlieren  study  of  the 
detonation  structure  of  regular  and  irregular  mixtures.  In  an  irregular  mix¬ 
ture,  2  CH4-2O2-0.2Air,  he  concluded  from  polar  calculations  based  on  the 
observed  wave  angles  that  the  transverse  wave  configurations  were  usually 
unreactive  and  reported  that  unburnt  islands  of  gas  were  formed.  In  these 
irregular  methane  mixtures,  the  unburnt  islands  were  fragmented  by  the  fine 
scale  cellular  structure.  These  detonations  were  reported  to  be  self-sustaining. 
In  a  regular  mixture,  2H2-C>2-3Ar,  he  reported  reactive  transverse  waves  and 
no  pocket  formation  except  when  the  initial  pressure  was  lowered  sufficiently 
for  the  detonation  to  become  marginal.  In  these  marginal  cases,  he  concluded 
both  reactive  and  unreactive  transverse  wave  configurations  were  observed. 
He  showed  schematically  how  regular,  triangular-shaped,  unburnt  gas  islands 
could  be  formed  by  the  collision  of  unreactive  transverse  waves.  He  noted  that 
in  these  mixtures,  the  detonation  was  unstable  and  failing.  Subbotin’s  con¬ 
clusions  about  the  reaction  zone  structure,  while  very  interesting,  are  based 
on  (often  unpublished)  schlieren  images,  and  furthermore,  result  from  exper¬ 
iments  conducted  in  a  very  narrow  (4.7  mm)  channel  which  was  also  rather 
short  (1.4  m  long),  making  it  difficult  to  determine  the  stability  of  the  deto¬ 
nation. 

In  1982,  Edwards  made  schlieren  images  in  a  2H2-02-4.5%Ar  mixture  in  a 
narrow  (6.3  mm)  channel  and  observed  density  gradients  of  irregular  regions 
of  gas  on  the  order  of  one  cell  width  downstream  of  the  front  [35].  An  ac¬ 
companying  numerical  simulation  of  a  detonation  near  its  initiation  limit  in 
this  mixture  shows  a  region  of  increased  induction  time  left  isolated  after  the 
collision  of  an  unreactive  triple  point  with  an  axis  of  symmetry.  In  more  re¬ 
cent  numerical  simulations,  Gamezo  et  al.  [36]  report  that  unburnt  pockets 
are  commonly  observed  and  conclude  mixtures  with  low  activation  energies 
(Ea/ RTvn— 2.1)  show  unburnt  pockets  that  are  triangular  in  shape  with  in¬ 
distinct  boundaries  while  higher  activation  energy  (Ea/ RTvn—7  A)  mixtures 
have  distinct,  irregular  pockets  which  survive  further  downstream.  Gamezo 
et  ah  [33]  believe  the  mechanism  by  which  the  pockets  are  consumed  also 
varies  with  the  activation  energy  of  the  mixture.  Pockets  in  low  activation  en¬ 
ergy  mixtures  appear  to  react  by  auto-ignition  after  shock  compression,  while 
pockets  in  irregular  mixtures  are  supposed  to  be  consumed  by  heat  and  mass 
exchange  with  neighboring  hot  gases  through  diffusion  and  turbulence. 

Sharpe  [37]  observes  the  collision  of  two  transverse  wave  configurations  in 
which  the  secondary  triple  point  becomes  detached  from  the  front  before  col¬ 
lision.  The  triple  point  configuration  that  remains  attached  to  the  front  is 
not  strong  enough  to  consume  the  unreacted  gas  upon  collision,  resulting  in 
an  unburnt  pocket.  Sharpe  also  shows  that  numerical  resolution  can  play  an 


2  stoichiometric  CH4-O2  diluted  with  6%  air 
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important  role  in  determining  the  structure  of  the  triple  point  and  the  pres¬ 
ence  of  unburnt  regions  with  insufficient  resolution  resulting  in  an  artificially 
accelerated  reaction  rate. 


All  previous  investigations  of  unreacted  pockets,  both  experimental  and  nu¬ 
merical,  have  been  carried  out  in  the  special  case  of  essentially  two-dimensional 
flows,  commonly  on  marginal  detonations.  Unburnt  pockets  have  been  postu¬ 
lated  to  appear  in  cases  where  the  transverse  waves  are  weak  and  do  not 
consume  the  reactants  behind  the  incident  wave  during  the  triple  point  col¬ 
lision.  Although  from  shock  polar  analysis  as  well  as  from  the  shape  of  the 
keystone  region  (see  section  below),  we  conclude  that  in  the  case  of  the  2H2- 
02-Ar  mixtures  we  studied,  the  transverse  waves  are  weak  (Mr  ~1.2)  and 
unreactive,  and  no  such  pockets  were  observed  in  the  present  experiments. 
The  PL1F  signal  intensity  decays  with  distance  from  the  front,  so  that  small 
differences  in  reactivity  may  be  difficult  to  image  more  than  2-3  cm  down¬ 
stream  of  the  reaction  zone.  However,  if  unreacted  pockets  exist,  we  would 
expect  to  see  evidence  of  their  formation  close  to  the  front.  The  PL1F  image 
from  shot  1417,  shown  in  detail  in  Fig.  14,  shows  what  appears  to  be  a  new 
Mach  stem  forming  just  after  a  triple  point  collision,  and  no  unreacted  region 
is  visible. 


In  the  N2-diluted  mixtures  studied,  some  unusual  features  are  observed  (see, 
for  example,  Shot  1427);  however,  these  do  not  have  the  structure  previously 
associated  with  unburnt  pockets.  We  attribute  these  features  to  the  relative 
orientation  of  the  light  sheet  to  the  planes  of  the  detonation  cells.  The  tran¬ 
sition  that  the  detonation  must  make  from  a  round  tube  to  the  square  test 
section  may  create  this  three-dimensional  effect.  A  study  with  sooted  foils  lo¬ 
cated  along  the  length  of  the  test  section  recorded  the  transition  in  cellular 
structure  from  spinning  to  orthogonal,  and  found  N2-diluted  mixtures  require 
a  greater  distance  for  the  cellular  structure  to  become  orthogonal,  resulting  in 
more  pronounced  three-dimensional  effects  in  these  mixtures. 


Sharpe’s  observations  about  numerical  resolution  [37]  may  explain  the  absence 
of  pockets  in  our  numerical  simulations.  However,  we  can  examine  the  poten¬ 
tial  for  pocket  formation  in  the  simulations  of  Eckett  [30]  by  computing  the 
shock  decay  time  through  numerical  differentiation  of  the  data  in  Fig.  12b. 
This  is  compared  to  the  computed  critical  decay  rates,  as  defined  by  Eckett  et 
al.  [34],  in  Fig.  13.  We  would  not  expect  any  decoupling  or  pocket  formation 
since  the  shock  decay  time  is  less  than  the  critical  shock  decay  time  through¬ 
out  the  entire  cell.  Also,  on  this  basis,  subsequent  computations  of  the  reaction 
zone  structure  will  use  the  quasi-steady  approximation  for  simplicity. 
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6  Triple  Point  Analysis 


Triple  point  configurations  through  a  detonation  cell  have  been  investigated 
experimentally,  theoretically,  and  numerically  by  various  researchers.  Progress 
has  been  made  towards  determining  the  triple  point  structure  at  different 
locations  through  the  cell  for  varying  mixtures  in  varying  geometries,  giving 
a  wide  spectrum  of  results. 

Because  of  the  difficulties  associated  with  studying  the  full  three-dimensional 
structure,  most  work  has  concentrated  on  marginal  detonation  waves  where 
the  detonation  is  propagating  close  to  its  limit,  usually  in  a  facility  in  which 
at  least  one  transverse  dimension  is  on  the  order  of  the  detonation  cell  width, 
A.  Experiments  have  concentrated  on  single-spin  detonation  [4,38,39]  or  on 
detonation  waves  in  a  narrow  channel  [4-8].  Two  general  classes  of  triple  point 
structure  have  been  suggested:  “weak”  configurations  that  involve  a  single 
triple  point  with  nonreactive  waves,  and  “strong”  conhgurations  that  involve 
a  second  triple  point,  generally  with  some  portion  of  the  transverse  wave 
becoming  reactive  (Fig.  15).  Many  variations  of  these  conhgurations  have  been 
reported.  As  discussed  in  Strchlow  and  Crooker  [7] ,  marginal  detonations  differ 
from  ordinary  detonations  in  several  ways;  in  particular,  the  transverse  wave 
spacing  and  strength  are  increased  in  comparison  to  ordinary  detonations.  For 
non-marginal  (or  ordinary)  detonations  in  three  dimensions,  the  triple  point 
structure  has  been  studied  only  indirectly,  predominantly  by  soot  foil  tracks 
or  with  visualization  techniques  that  integrate  through  the  how  held. 

The  strength  and  reactivity  of  transverse  waves,  at  least  in  the  vicinity  of 
the  triple  point,  may  be  calculated  by  polar  analysis  and  zero-dimensional 
reaction  zone  models,  as  discussed  in  detail  later  in  this  section.  From  this 
type  of  analysis,  we  can  order  some  sample  mixtures  in  terms  of  their  expected 
transverse  wave  strengths  so  as  to  have  a  loose  classification  of  the  mixtures 
studied  in  the  present  experiments. 

It  should  be  noted  that  the  calculations  presented  in  Table  2  are  all  for  a  con¬ 
stant  track  angle  (0=33°),  and,  as  will  be  discussed  later,  that  the  transverse 
wave  strength  is  a  strong  function  of  the  track  angle.  However,  in  a  study  of 
the  experimentally  obtained  track  angle  for  different  mixtures,  Strchlow  and 
Biller  [40]  found  a  marked  variation  in  the  track  angles  present  on  a  single 
foil  but  no  detectable  relationship  between  track  angle  and  initial  pressure 
mixture  or  dilution.  Their  study  contributed  to  our  decision  to  keep  the  track 
angle  constant  in  our  calculations. 

In  general,  we  calculate  that  Ar-diluted  mixtures  have  weaker  transverse  waves 
than  N2-diluted  mixtures  which  are  weaker  in  turn  than  mixtures  with  C02 
dilution  for  the  same  percent  diluent.  Hydrocarbons  have  stronger  transverse 
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waves  than  H2,  with  transverse  wave  strength  increasing  with  carbon  content 
for  methane,  ethylene,  and  propane.  Strehlow  [41]  made  a  soot  foil  study 
of  cellular  regularity  for  various  mixtures,  and  we  find  that  the  calculated 
transverse  wave  strength,  (P3-P2) / P2,  can  be  related  to  the  cellular  regularity. 
Strehlow  observed  very  regular  cellular  structure  in  Ar-diluted  mixtures.  The 
cellular  regularity  of  2H2-O2  mixtures  increased  with  increasing  Ar  dilution, 
while  there  was  no  apparent  effect  up  to  50%  N2  dilution.  The  addition  of 
C02  was  found  not  to  improve  the  regularity.  Ethylene-oxygen  soot  foils  were 
reported  to  be  considerably  more  complex  than  hydrogen-oxygen  foils.  The 
structural  regularity  of  CH4-2O2  and  C3H8-502  mixtures  was  classified  as 
poor.  Although  the  variation  in  the  calculated  transverse  wave  strength  for 
these  mixtures  is  small  (the  normalized  pressure  jump  varies  from  0.33  to  0.40), 
the  variation  in  cellular  regularity  observed  in  soot  foil  tracks  is  noticeable. 
Comparing  Strehlow’s  results  with  the  calculations  in  Table  2,  it  can  be  seen 
that  cellular  regularity  decreases  with  increasing  transverse  wave  strength. 

The  reactivity  of  the  transverse  waves  can  be  ordered  by  77/77,  the  ratio  of 
the  constant  volume  induction  time  for  a  particle  passing  through  the  incident 
wave  and  then  immediately  through  the  transverse  wave,  77,  normalized  by 
the  induction  time  behind  the  incident  wave,  77.  This  parameter  has  a  range 
0  <  77/77  <  1.  Values  close  to  0  correspond  to  a  very  short  induction  time  be¬ 
hind  the  transverse  wave,  so  that  the  transverse  wave  may  be  called  reactive. 
Values  close  to  1  correspond  to  a  relatively  long  induction  time,  meaning  the 
transverse  wave  has  a  negligible  effect  on  the  reaction.  For  all  the  mixtures 
shown  here,  77/77  is  about  0.5.  Our  calculations  of  the  induction  time  behind 
the  transverse  wave,  77,  show  a  short  induction  time  in  the  case  of  2H2-O2- 
3Ar  and  a  longer  time  in  CH4-2O2.  This  trend  agrees  with  the  observations 
of  Subbotin  [8],  who  classified  the  transverse  waves  in  2H2-02-3Ar  mixtures 
as  reactive  and  the  transverse  waves  in  CH4-2O2-0.2Air  as  unreactive.  As  can 
be  seen  from  the  table,  the  mixtures  considered  in  this  study  have  transverse 
waves  with  a  Mach  number  of  about  1.3  and  strength,  defined  as  the  normal¬ 
ized  pressure  jump  across  the  wave,  (P3  —  P2)/TV  of  about  0.35.  (Since  for  an 
oblique  wave  only  the  normal  component  of  the  flow  is  affected,  the  normalized 
pressure  jump  is  a  better  measure  of  the  strength  of  the  transverse  waves  than 
the  My.)  Strehlow  and  Biller  [40]  found  wave  strengths  ranging  from  0.22  to 
1.2  were  possible,  so  the  transverse  waves  in  the  present  experiments  are  at 
the  lower  end  of  this  range. 

Shock  and  detonation  polar  calculations  are  carried  out  to  analyze  the  triple 
point  configurations  observed  in  the  PLIF  images.  This  technique  has  been 
used  previously  by  several  researchers  [42-44,8].  The  shock  polar  is  calculated 
using  the  oblique  shock  jump  relations  and  assumes  a  perfect  gas,  with  the 
ratio  of  specific  heats  taken  to  be  that  of  the  reactants,  74.  The  detonation 
polar  [45]  is  a  locus  of  solutions  possible  for  an  oblique  detonation  which  is 
overdriven  at  shock  angles  greater  than  that  corresponding  to  the  CJ  wave.  A 
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two-gamma  model  is  used  to  model  the  variation  in  thermodynamic  properties 
between  the  reactants  and  products.  The  analysis  assumes  the  flow  is  quasi- 
steady  in  the  frame  of  the  triple  point  and  the  waves  are  straight  in  the  vicinity 
of  the  triple  point,  see  Fig.  16.  This  assumption  is  more  valid  for  a  nonreactive 
flow  than  reactive,  and  a  probable  effect  of  the  chemistry  will  be  to  produce 
some  wave  curvature  in  the  vicinity  of  the  triple  point. 

The  incident  shock  Mach  number  at  each  location  was  estimated  from  the  cal¬ 
culated  centerline  pressure  profile  discussed  in  the  previous  section  (Fig.  12a) 
and  the  incident  flow  angle  was  taken  to  be  57°  (corresponding  to  a  track 
angle,  0,  of  33°).  This  value  is  representative  of  track  angles  recorded  on  soot 
foils  mounted  on  the  side  walls  of  the  test  section.  Given  the  incident  wave 
Mach  number  and  the  track  angle,  the  triple  point  configuration  can  be  solved 
by  matching  the  pressure  and  flow  deflection  angle,  6,  across  the  slip  stream. 
The  shock  and  detonation  polars  are  therefore  conveniently  displayed  on  a 
pressure-deflection  (P-0)  plot.  An  example  of  such  a  calculation  together  with 
the  resulting  triple  point  configuration  for  a  stoichiometric  H2-O2  mixture  with 
70%  Ar  dilution  at  6.67  kPa  initial  pressure  is  shown  in  Fig.  17.  Note  that  as 
there  is  no  intersection  of  the  detonation  and  transverse  wave  polars,  there  is 
no  solution  possible  which  involves  a  reactive  Mach  stem,  i.e.  the  existence  of 
an  induction  zone,  however  small,  is  crucial  for  obtaining  a  solution.  For  this 
and  all  other  mixtures  considered  in  this  study,  a  solution  consisting  of  a  single 
triple  point  which  involved  the  detonation  polar  for  the  leading  wave  was  not 
possible;  thus,  the  Mach  stem  is  always  modeled  as  non-reactive.  In  a  study  of 
an  artificially-produced  reactive  Mach  stem  from  interferograms  obtained  by 
White  and  Caby  [46],  Strchlow  [47]  concludes  that  a  non- reactive  triple  point 
balance  is  appropriate  in  the  region  of  the  primary  triple  point. 

A  single  non-reactive  slip  stream  balance  for  the  2H2-C>2-7Ar  mixture  at 
6.67  kPa  is  compared  to  numerical  schlieren  images  (Fig.  18)  and  found  to 
be  in  reasonable  agreement  (to  within  the  resolution  of  the  calculation)  close 
to  the  triple  point.  As  can  be  seen  from  the  figure,  the  greatest  discrepancy  is 
in  the  transverse  wave  angle  away  from  the  triple  point;  the  transverse  wave  is 
calculated  to  be  almost  perpendicular  to  the  incident  wave  by  the  local  non¬ 
reactive  analysis.  However,  the  transverse  wave  will  interact  with  the  reaction 
zone  behind  the  incident  wave,  curving  the  wave  towards  the  cell  centerline  as 
it  moves  into  a  region  with  higher  speed  of  sound.  An  alternative  explanation 
is  the  existence  of  a  second  triple  point,  unresolved  by  the  numerical  simula¬ 
tion.  The  single  triple  point  calculation  does  not  preclude  the  existence  of  a 
more  complex  Mach  reflection  since  a  second  triple  point  calculation  does  not 
affect  the  initial  triple  point  solution,  but  the  angle  of  the  downstream  portion 
of  the  transverse  wave  is  likely  to  be  affected  by  such  a  secondary  triple  point. 

An  attempt  was  made  to  use  shock  and  detonation  polars  to  calculate  a  triple 
point  configuration  involving  secondary  triple  points.  Following  the  procedure 
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commonly  used  to  calculate  double  Mach  reflections,  the  calculation  was  first 
made  in  the  frame  of  the  first  triple  point  as  discussed  previously,  then  a 
transformation  was  made  to  the  coordinate  system  of  the  second  triple  point 
where  a  similar  second  triple  point  calculation  can  be  made.  To  perform  this 
shift,  the  Law-Glass  [48]  model  was  used.  This  model  assumes  the  second  triple 
point  moves  with  the  velocity  induced  in  the  flow  behind  the  incident  wave 
(region  2,  see  Fig.  16)  in  the  laboratory  frame.  For  all  mixtures  considered 
in  the  study,  the  flow  in  region  3  becomes  subsonic  under  this  transformation 
and  a  double  triple  point  configuration  is  not  possible.  However,  the  Law- 
Glass  model  was  conceived  for  unreactive  flow  over  a  wedge  and  its  use  here 
is  certainly  speculative,  so  the  existence  of  a  secondary  triple  point  remains 
an  open  question. 

Single  triple  point  configurations  were  calculated  for  the  mixtures  in  the  ex¬ 
periments:  stoichiometric  H2-O2  with  80%  Ar,  with  85%  Ar,  and  with  60%  N2 
dilution,  all  at  20  kPa  initial  pressure.  The  calculated  values  for  stoichiometric 
H2-O2  with  70%  Ar  at  6.67  kPa  is  also  shown  for  comparison.  Results  for  one 
location  in  the  cell  are  shown  in  Table  3. 

The  Mach  stem  velocity,  Mmh,  shown  above  does  not  correspond  to  the  calcu¬ 
lated  centerline  Mach  stem  velocity  (Fig.  12b),  which  implies  the  Mach  stem 
is  curved.  Except  immediately  at  the  apex  of  the  cell,  the  centerline  velocity 
of  the  Mach  stem  is  lower  than  that  predicted  by  the  local  triple  point  anal¬ 
ysis,  indicating  the  wave  is  weaker  away  from  the  triple  point.  The  induction 
distance  is  therefore  greatest  at  the  centerline  and  the  Mach  stem  shock  may 
have  more  curvature  than  is  apparent  in  the  PLIF  images  of  the  OH  contour. 

The  transverse  wave  strength  and  shock  angle  are  very  sensitive  to  changes 
in  the  track  angle.  From  sooted  foils  mounted  on  the  side  wall  of  the  test 
section,  we  observe  ccll-to-cell  variations  in  the  track  angle  of  around  ±3°  for 
each  mixture.  The  ccll-to-cell  variation  results  in  a  range  of  transverse  wave 
strengths.  In  order  to  make  qualitative  comparisons  of  the  triple  point  con¬ 
figurations  in  different  mixtures,  the  track  angle  was  held  constant  for  the 
analysis  shown  above.  A  study  of  the  influence  of  track  angle  is  presented  in 
Table  4  where  the  minimum,  maximum,  and  average  track  angles  measured 
from  the  side  wall  soot  foils  are  used  in  calculating  the  triple  point  configura¬ 
tion.  As  can  be  seen  from  the  table,  changes  in  track  angle  have  a  significant 
effect  on  the  transverse  wave  structure.  It  can  also  be  noted  that  since  the 
incident  wave  velocity  does  not  change  a  great  deal  through  the  second  half  of 
the  cell  (Fig.  12b),  changes  in  the  strength  and  shock  angle  of  the  transverse 
wave  through  the  cell  may  be  attributed  to  variations  in  the  track  angle.  An 
increase  in  track  angle  through  the  second  half  of  the  cell  is  evident  in  soot 
foil  tracks  of  some  detonations,  particularly  those  with  orthogonal  structure 
[41],  As  the  transverse  wave  strength  increases  with  increasing  track  angle, 
this  is  consistent  with  the  evolution  of  the  transverse  wave  structure  observed 
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by  some  authors  [7,49]. 


In  general,  more  detailed  information  on  all  of  these  quantities:  incident  and 
Mach  wave  strength,  triple  point  track  angle,  and  wave  curvature  at  the  triple 
point  will  be  required  before  sharper  bounds  on  the  triple  point  configuration 
can  be  obtained.  Experimentally,  this  has  been  difficult  to  do  because  of  the 
three-dimensional  nature  of  the  cellular  structure.  Numerically,  few  computa¬ 
tions  have  been  carried  out  with  the  spatial  resolution  or  sufficiently  realistic 
reaction  mechanism  in  order  to  confidently  obtain  these  properties. 


Constant  volume  explosion  calculations  were  made  to  obtain  OH  mole  fraction 
contours  in  the  region  of  the  triple  point.  A  detailed  chemical  mechanism 
[50],  previously  validated  against  shock  tube  ignition  delay  data  [28],  was 
used.  The  constant  volume  calculations  were  carried  out  for  flow  along  several 
particle  paths  through  the  calculated  triple  point  configuration  (see  Fig.  19).  A 
contour  is  drawn  linking  the  locations  along  each  path  at  which  the  OH  mole 
fraction  reaches  a  specified  value  representative  of  the  end  of  the  induction 
zone.  Distance  along  the  path  is  calculated  by  multiplying  the  time  by  the 
post-shock  velocity.  Particle  path  1  passes  directly  through  the  Mach  stem. 
The  flow  along  particle  paths  2-10  passes  through  the  incident  wave  and  is 
partially  reacted  before  passing  through  the  transverse  wave.  The  transverse 
wave  was  assumed  to  remain  straight  with  constant  Mach  number.  Particle 
path  11  passes  through  only  the  incident  wave  before  the  OH  mole  fraction 
reaches  the  defined  value.  The  dependance  of  the  location  of  the  rise  in  OH 
mole  fraction,  and  therefore  of  the  induction  distance,  on  the  varying  shock 
strength  in  different  parts  of  the  front,  results  in  the  discontinuous  keystone 
structure  we  see  in  the  PLIF  images. 


Fig.  19a  shows  the  calculated  OH  mole  fraction  contour  for  a  transverse  wave 
strength  of  0.325  in  2H2-02-17Ar,  initial  pressure  20  kPa.  Other  mixtures 
considered  in  experiments  have  similar  profiles.  Varying  the  track  angle  has 
the  most  significant  effect  on  the  OH  contour.  Fig.  19b  shows  a  calculated  OH 
mole  fraction  contour  in  the  same  mixture  for  a  significantly  increased  track 
angle,  0  =  40°,  corresponding  to  a  transverse  wave  strength  of  0.666. 


These  calculations  show  quantitative  agreement  with  the  PLIF  images,  al¬ 
though  the  location  of  the  shock  structure  must  be  estimated,  and  with  nu¬ 
merical  simulations.  Results  show  the  transverse  waves  do  not  significantly 
accelerate  the  reaction,  and  the  shear  layers  separate  reacted  and  unreacted 
gas  and  form  the  keystones. 
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6. 1  The  role  of  transverse  waves 


A  key  scientific  question  is:  What  is  the  role  of  the  transverse  waves  in  deto¬ 
nation  propagation?  Two  extreme  points  of  view  are: 

(1)  The  transverse  waves  are  simply  manifestations  of  instability  and  do  not 
play  an  essential  role  in  the  propagation  mechanism.  The  main  portion  of 
the  reaction  is  due  to  the  leading  shock  waves  and  only  a  small  fraction 
of  the  reactants  pass  through  the  transverse  wave. 

(2)  Transverse  waves  are  essential  to  detonation  propagation.  The  transverse 
waves  and  interactions  of  transverse  waves  provide  high-temperature  re¬ 
gions  (“hot  spots”)  that  serve  as  reaction  centers.  A  substantial  amount 
of  the  chemical  reaction  takes  place  behind  the  transverse  waves  or  is 
associated  with  reaction  centers. 

We  have  examined  only  one  fuel-oxidizer  combination  and  two  diluents  in  the 
present  study.  For  those  mixtures,  we  observe  behavior  that  is  closer  to  the  first 
point  of  view:  the  transverse  waves  play  a  nonessential  role  in  detonation  prop¬ 
agation.  Mixtures  with  very  different  thermodynamic  and  kinetic  parameters 
need  to  be  examined  to  determine  if  the  second  type  of  behavior  is  possible. 
On  the  basis  of  experiments  in  porous  wall  tubes,  Radulescu  and  Lee  [51]  have 
proposed  that  mixtures  with  high  activation  energies,  associated  with  certain 
hydrocarbon  fuels  without  inert  gas  dilution,  are  more  representative  of  the 
second  class  of  behavior.  Further  studies  are  in  progress  [52]  in  onr  laboratory 
to  examine  a  wider  range  of  mixtures  and  determine  the  chemical  structure 
of  the  detonation  front  for  those  cases. 

The  present  study  was  restricted  to  mixtures  with  modest  activation  energy 
(see  Table  1)  but  included  mixtures  with  a  range  of  cellular  regularity.  The 
highly  argon-diluted  mixtures  are  well  known  to  have  “excellent  regularity” 
while  the  nitrogen-diluted  mixtures  are  classified  [53]  as  having  “fair”  regu¬ 
larity  (see  Fig.  9).  ft  is  important  to  note  that  our  OH  PL1F  visualizations  of 
these  two  mixtures  do  not  show  a  dramatic  difference  in  the  structure  of  the 
front  that  would  indicate  a  qualitative  difference  in  the  role  of  the  transverse 
waves.  We  do  see  a  more  irregular  structure  of  the  OH  front  in  the  nitrogen- 
diluted  cases  than  in  the  argon-diluted  cases,  but  we  do  not  see  any  indication 
of  pockets  or  inclusions  of  unburnt  material  behind  the  leading  shock  front. 
We  conclude  that  the  irregularity  is  simply  due  to  the  spontaneous  generation 
and  decay  of  triple  points  along  the  front,  and  the  transverse  waves  play  a 
similar  role  for  both  diluents. 

From  onr  simple  considerations  about  triple  point  structure,  we  find  a  mod¬ 
est  decrease  in  the  reaction  time  behind  the  transverse  wave  structures  that 
are  consistent  with  onr  present  observations.  The  reaction  time  for  the  fluid 
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elements  closest  to  the  triple  point  will  decrease  about  a  factor  of  two  due 
to  the  increase  in  temperature  across  the  transverse  wave.  While  this  affects 
the  details  of  the  keystone  shape,  it  does  not  seem  appropriate  to  refer  to 
the  transverse  waves  as  “reactive”  in  this  case.  However,  it  is  clear  that  the 
possibility  of  reactive  transverse  waves  has  been  demonstrated  by  computa¬ 
tions  and  experiments  in  marginal  propagation  situations.  Further  studies  are 
needed  with  other  mixtures  to  determine  if  reactive  transverse  waves  have  a 
role  in  “ordinary”  (near  CJ)  detonations. 


7  Conclusions 


Experimental  observations  have  been  made  of  the  reaction  zone  structure  in 
propagating  detonations.  H2-02-diluent  mixtures  representative  of  regular  and 
irregular  cell  structures  were  studied.  PLIF  images  reveal  distinctive  keystone 
features  associated  with  the  triple  points.  Differences  in  cellular  structure  reg¬ 
ularity  are  seen  in  the  different  mixtures,  with  N2-diluted  mixtures  sometimes 
exhibiting  features  that  may  be  attributed  to  three-dimensional  effects. 

The  local  structure  of  the  triple  point  is  analysed  by  an  unreactive  shock 
polar  calculation,  revealing  that  keystone  features  are  due  to  a  discontinuity 
in  reactivity  across  the  shear  layer.  From  this  analysis  and  from  the  images 
themselves,  we  conclude  that  transverse  waves  have  little  influence  on  the 
reaction  in  the  case  of  the  Ar-diluted  mixtures. 

Previous  numerical  work  has  suggested  that  in  the  case  of  mixtures  with  low- 
amplitude  transverse  waves,  pockets  of  unreacted  gas  may  be  formed  if  the 
transverse  waves  do  not  consume  the  unreacted  gas  behind  the  incident  wave 
during  triple  point  collision.  A  region  of  unreacted  gas  becomes  isolated  from 
the  front  and  may  persist  a  cell  length  downstream.  In  all  images  obtained,  no 
detached  pockets  of  gas  are  seen  within  the  held  of  view  of  the  measurement. 

Our  study  has  examined  a  limited  range  of  mixtures  and  initial  conditions. 
Other  triple  point  configurations,  reactive  transverse  waves,  and  unreacted 
pockets  may  be  possible  for  different  mixtures  with  a  wider  range  of  chemical 
and  thermodynamic  properties.  Experiments  are  in  progress  to  examine  these 
situations. 
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Table  Captions 


Table  1:  Mixture  parameters  for  the  mixtures  considered  in  this  study,  calcu¬ 
lated  with  STANJAN  [27].  P\  is  the  initial  pressure.  71  is  the  ratio  of  specific 
heats  of  the  reactants  that  was  used  in  calculating  the  shock  polars  with  the 
perfect  gas  assumption.  The  induction  zone  length,  ACT,  is  defined  by  the  point 
of  maximum  heat  release;  Ea/RTV^  is  the  non-dimensional  activation  energy, 
calculated  as  in  [28].  Both  these  parameters  are  calculated  using  a  detailed 
chemical  mechanism  [29].  “:ccll  widths  measured  in  test  section. 

Table  2:  Mixture  parameters  and  transverse  wave  strengths  for  some  sample 
mixtures,  calculated  in  the  vicinity  of  a  primary  triple  point.  All  calculations 
are  for  mixtures  at  20  kPa.  The  incident  wave  Mach  number  is  kept  constant  at 
Mcj-  The  transverse  wave  strength  is  not  a  strong  function  of  initial  pressure 
until  the  detonation  becomes  marginal.  Activation  energy  and  induction  time 
calculations  use  the  detailed  chemical  mechanism  of  Konnov  [29],  previously 
validated  against  shock  tube  ignition  delay  data  [28]. 

Table  3:  Triple  point  configuration  at  x  —  xa=  0.7 L,  ( U /Ucj=0.94: ),  with 
0  =  33°. 

Table  4:  Triple  point  configuration  with  varying  track  angle  for  2H2-02-12Ar 
at  20  kPa  initial  pressure.  Triple  point  values  are  the  minimum,  maximum 
and  average  value  observed  on  a  soot  foil  for  this  mixture.  The  incident  wave 
Mach  number  is  kept  constant  at  Mcj- 
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Table  1 


Composition  Pi  71  Ucj  Tvn 


(kPa) 

(m/s) 

(K) 

2H2-02-7Ar 

6.67 

1.568 

1618 

1904 

2H2-02-12Ar 

20 

1.602 

1518 

1899 

2H2-02-17Ar 

20 

1.622 

1415 

1775 

2H2-O2-20Ar 

20 

1.630 

1357 

1687 

2H2-02-4.5N2 

20 

1.405 

1874 

1454 

2H2-02-5.6N2 

20 

1.405 

1797 

1403 

Table  2 

Mixture 

C'PvN 

kJ /kg-K 

Cl 

m/s 

MT  (Pi  -  P2)/P2 

Ea/RTvN 

tt 

Us 

tt/ti 

2H2-02 

2.87 

537 

1.4 

0.38 

5.4 

0.25 

0.54 

2H2-02-3Ar 

1.08 

379 

1.3 

0.36 

5.0 

0.32 

0.54 

2H2-02-12Ar 

0.69 

339 

1.2 

0.33 

5.2 

0.94 

0.51 

2H2-02-3N2 

1.72 

416 

1.4 

0.38 

6.1 

0.89 

0.52 

2H2-02-C02 

1.96 

411 

1.4 

0.38 

10.1 

1.5 

0.54 

CH4-202 

2.09 

355 

1.5 

0.38 

11.8 

1.3 

0.44 

C2H4-302 

1.86 

326 

1.6 

0.39 

6.8 

0.38 

0.58 

C3H8-502 

2.01 

306 

1.7 

0.40 

10.5 

0.28 

0.56 

Figure  Captions 


Figure  1:  a)  Schematic  of  the  detonation  tube  with  attached  “cookie-cutter” 
and  test  section,  b)  Schematic  of  the  experimental  set  up  for  simultaneous 
use  of  schlieren  and  PLIF  techniques.  In  experiments  where  only  PLIF  im¬ 
ages  were  acquired,  the  camera  is  placed  exactly  perpendicular  to  the  light 
sheet.  The  flow  transitions  from  the  round  tube  to  the  square  test  section 
with  optical  access.  In  experiments  without  flow  visualization,  soot  foils  can 
be  mounted  both  on  the  square  end  wall  of  the  test  section  and  on  the  side 
wall. 

Figure  2:  Images  of  OH  fluorescence  behind  detonation  fronts  in  stoichiomet¬ 
ric  hydrogen-oxygen  mixture  diluted  with  80%  argon,  initial  pressure  20  kPa. 

Figure  3:  Images  of  OH  fluorescence  behind  detonation  fronts  in  stoichiomet¬ 
ric  hydrogen-oxygen  mixture  diluted  with  85%  argon,  initial  pressure  20  kPa. 

Figure  4:  Images  of  OH  fluorescence  behind  detonation  fronts  in  stoichiometric 
hydrogen-oxygen  mixture  diluted  with  60%  (shot  1426)  or  65%  (shots  1427 
and  1428)  nitrogen,  initial  pressure  20  kPa. 

Figure  5:  OH  number  density  and  temperature  ZND  profiles  for  a  CJ  detona¬ 
tion  in  2H2-02-17Ar  at  20  kPa  initial  pressure.  The  solid  line  corresponds  to 
OH  number  density,  the  dashed  line  to  the  temperature  profile. 

Figure  6:  Two  soot  foils  (150x150  mm)  obtained  from  the  end  wall  of  the 
test  section  in  separate  shots,  both  are  2H2-02-12Ar  at  20  kPa.  Note  that  the 
transverse  wave  families  are  not  necessarily  orthogonal  to  each  other  or  to  the 
test  section  walls. 

Figure  7:  Schematic  illustration  of  possible  three-dimensional  effects  due  to 
non-orthogonal  intersection  of  the  transverse  wave  systems  and  the  laser  light 
sheet. 

Figure  8:  Simultaneous  schlieren  and  OH  fluorescence  images  behind  detona¬ 
tion  fronts  in  stoichiometric  hydrogen-oxygen  mixture  diluted  with  85%  argon 
(shots  1432  and  1433)  and  87%  argon  (shot  1434),  initial  pressure  20  kPa.  (a) 
Overlay  of  PLIF  and  schlieren  images  (b)  schlieren  (c)  PLIF. 

Figure  9:  Soot  foils  obtained  in  a  280  mm  diameter  detonation  tube  for  2H2- 
02-12Ar  (a)  and  2H2-02-4.5N2  (b),  both  at  20  kPa  initial  pressure.  Flow  is  top 
to  bottom.  The  scale  shown  corresponds  to  50  mm.  The  Ar-diluted  mixture 


is  an  example  of  regular  cell  structure,  while  the  N2-diluted  mixture  has  more 
irregular  cell  structure. 

Figure  10:  Results  of  two-dimensional  numerical  simulation  of  Chapman-Jou- 
guet  detonation  in  a  60  mm  wide  channel:  a)  Numerical  schlieren  b)  OH  mass 
fraction  distribution.  The  mixture  is  stoichiometric  hydrogen-oxygen  diluted 
with  70%  argon,  initial  pressure  6.67  kPa. 

Figure  11:  (a)  Close-up  of  OH  mass  fraction  contours  shown  in  Fig.  10.  (b) 
Close  up  of  schlieren  image  of  density  field,  (c)  Schematic  of  detonation  front 
features. 

Figure  12:  (a)  Pressure  behind  the  leading  shock  wave  along  the  centerline  of  a 
detonation  cell  from  the  simulations  of  Fig.  10.  These  results  were  obtained  by 
interrogation  of  the  pressure  field  in  the  fully  developed  regime  (see  Fig.  3.39 
in  Eckett  [30]).  (b)  The  Mach  number  of  the  leading  shock  front  along  the 
centerline  of  the  detonation  cells  for  the  simulations  of  Fig.  10.  Velocities  were 
obtained  by  differentiation  of  the  smoothed  wave  trajectories  obtained  by  in¬ 
terrogating  the  solution  held  in  the  fully  developed  regime  (see  Fig.  3.37  in 
Eckett  [30].) 

Figure  13:  Time  scales  in  detonation  cell  for  simulations  of  Fig.  10.  The  shock 
decay  times  were  obtained  by  numerical  differentiation  of  smooth  data  shown 
in  Fig.  12b.  The  induction  times  are  calculated  from  ZND  simulations  based 
on  the  instantaneous  values  of  the  leading  shock  velocity.  The  critical  decay 
time  concept  is  discussed  in  Eckett  et  al.  [34], 

Figure  14:  Detail  from  shot  1417,  2H2-02-17Ar  at  20  kPa,  showing  a  new  Mach 
stem  formed  after  a  triple  point  collision.  No  unreacted  pockets  are  evident. 

Figure  15:  Cartoon  to  illustrate  the  nomenclature  associated  with  the  triple 
point.  Two  triple  points  are  shown:  a  primary  triple  point  at  the  intersection 
of  the  Mach  stem  and  incident  wave,  and  a  secondary  triple  point  downstream 
of  the  main  front.  The  how  downstream  of  the  second  triple  point  is  shown 
here  as  unreactive,  but  some  authors  describe  a  reaction  zone  closely  coupled 
to  the  portion  of  the  transverse  wave  downstream  of  the  second  triple  point, 
the  so-called  “strong”  triple  point  configuration.  A  “weak”  configuration  gen¬ 
erally  omits  the  second  triple  point. 

Figure  16:  Wave  angles  for  shock  polar  calculation. 

Figure  17:  Shock  and  detonation  polar  calculations  for  stoichiometric  hydrogen- 
oxygen  mixture  diluted  with  70%  argon,  initial  pressure  6.67  kPa,  x  —  xD= 
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0.91L,  0=33°. 


Figure  18:  Triple  point  polar  calculation  from  Fig.  17  superimposed  on  numeri¬ 
cal  schlieren  for  x—x0=  0.91L,  in  stoichiometric  hydrogen-oxygen  diluted  with 
70%  argon,  initial  pressure  6.67  kPa. 

Figure  19:  Prediction  of  the  keystone  region  shape  for  stoichiometric  hydrogen- 
oxygen  with  85%  argon  dilution  for  0  =  33°  (a)  and  0  =  40°  (b).  These 
predictions  are  based  on  the  idealized  model  of  the  triple  point  configuration 
discussed  in  the  text  and  estimates  of  the  OH  mole  fraction  using  simple  zero 
and  one-dimensional  reaction  zone  models.  The  transverse  wave  was  assumed 
to  remain  straight  with  constant  Mach  number.  The  edge  of  the  reaction  zone 
is  defined  as  the  location  of  the  contour  xoh=0.005.  The  choice  of  OH  contour 
Xo//=0.005  is  arbitrary,  but  as  the  OH  rise  is  exponential,  the  features  remain 
qualitatively  the  same  irrespective  of  which  value  is  chosen. 
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Fig.  1.  one  column  images 
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Fig.  2.  two  column  image 
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Fig.  3.  two  column  image 
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Fig.  4.  two  column  image 
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Fig.  5.  one  column  image 
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Fig.  7.  one  column  image 
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Fig.  8.  two  column  images 
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Fig.  9.  two  column  image 
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Fig.  10.  two  column  image 
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Fig.  11.  two  column  image 
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Fig.  13.  one  column  image 
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Fig.  14.  one  column  image 
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Fig.  15.  one  column  image 
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Fig.  16.  one  column  image 
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Fig.  17.  one  column  image 
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Fig.  18.  one  column  image 
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Fig.  19.  one  column  image 
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